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We discuss recent advances in the calculation of thermodynamic observables using improved actions. In particular, 
we discuss the calculation of the equation of state of the SU{3) gauge theory, the critical temperature in units of the 
string tension, the surface tension and the latent heat at the deconfinement transition. We also present first results 
from a calculation of the equation of state for four-flavour QCD using an O(a^) improved staggered fermion action 
and discuss possible further improvements of the staggered fermion action. 


1. Introduction 

The idea that strongly interacting hadronic mat¬ 
ter undergoes a phase transition to a new phase, 
the quark-gluon plasma, has been around for a long 
time. More than 10 years ago lattice calculations 
have given first direct evidence for the existence of 
such a phase transition and opened the way for its 
detailed quantitative analysis. Many of the non- 
perturbative features of finite temperature QCD, 
which originally emerged from perturbative stud¬ 
ies have since then been analyzed in lattice calcu¬ 
lations. However, many of these numerical studies 
also had to remain on a qualitative level for a long 
time because simulations on coarse lattices were 
hampered by large discretization errors and calcu¬ 
lations at smaller lattice spacing suffered from poor 
statistical accuracy. 

The discretization scheme proposed originally by 
K. Wilson for gauge theories does for a finite 
lattice cut-off, a, introduce systematic O(a^) devi¬ 
ations from the continuum formulation. It is well 
known that this introduces severe problems in ther¬ 
modynamic calculations, which in the past made a 
direct determination of the physics in the contin¬ 
uum limit difficult. Already the thermodynamics of 
free Bose (gluon) or Fermi (quark) gases deviates 
strongly from the continuum ideal gas results when 
calculated on coarse Euclidean lattices This 
problem carries over to QCD, where bulk thermo¬ 
dynamic observables like energy density and pres¬ 
sure do rapidly come close to the non-interacting 
ideal gas limit above the deconfinement phase tran¬ 
sition. The calculation of these quantities is thus 
expected to suffer from similarly strong cut-off ef¬ 
fects as the ideal gas, despite the fact that the high 
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temperature plasma phase does in many respects 
still show strong non-perturbative properties, char¬ 
acterized e.g. by thermal screening masses and 
quasi-particle excitations. 

At high temperature the relevant contributions 
to thermodynamic observables result from mo¬ 
menta which are of the order of the temperature, 
i.e. (p) ^ 3T. However, in lattice calculations T is 
fixed through the temporal extent. Nr, of the lat¬ 
tice and the cut-off, a, i.e. T = 1/Nra. As the com¬ 
putational effort in the calculation of bulk thermo¬ 
dynamic observables, which are of dimension 
increases approximately like , it is evident that 
the temporal extent of the lattice has to remain 
rather small in most thermodynamic calculations. 
The relevant momenta for e.g. lattices with tempo¬ 
ral extent Nr = 4 are therefore close to the cut-off 
where the lattice and continuum dispersion rela¬ 
tions differ strongly from each other. Indeed this is 
the origin of the well-known discrepancy between 
the energy density of an ideal gas calculated on a 
finite lattice {e{Nr)) and in the continuum (csb) 
Q. The cut-off dependence in finite temperature 
calculations thus shows up as a finite size effect, 
which should not be confused with the finite size 
dependence resulting from the spatial extent of the 
lattice. The latter is an infra-red effect and controls 
the approach to the thermodynamic limit. 

In the next section we briefly discuss some im¬ 
proved actions, which recently have been used 
in thermodynamic calculations. In Section 3 we 
present results from a perturbative calculation of 
cut-off effects in the ideal gas (infinite tempera¬ 
ture) limit. Results from numerical simulations of 
the SU{3) gauge theory and four-flavour QCD are 
discussed in Sections 4 and 5. In Section 6 we sug¬ 
gest a new staggered fermion action suitable for 
thermodynamic calculations. Finally we give our 
conclusions in Section 7. 
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2. Improved actions 

During the last few years much progress has been 
made in dealing with the systematic discretization 
errors in lattice regularized quantum field theories. 
Various improved discretization schemes for the La- 
grangian of QCD have been constructed and ex¬ 
plored, which do show much less cut-off dependence 
than the prescription originally given by Wilson 

§■ 

When formulating a discretized version of QCD 
one has a great deal of freedom in choosing a lattice 
action. Different formulations may differ by sub¬ 
leading powers of the lattice cut-off, which vanish 
in the continuum limit. This has, for instance, been 
used by Symanzik to systematically improve scalar 
field theories and has then been applied to lat¬ 
tice regularized SU{N) gauge theories |^,^. In ad¬ 
dition to the elementary plaquette term appearing 
in the standard Wilson formulation of lattice QCD 
larger loops can be added to the action in such a 
way that the leading 0{a^g^) deviations from the 
continuum formulation are eliminated and correc¬ 
tions only start in 0{a'^g°, a^g^). A simple class 
of improved actions is, for instance, obtained by 
adding planar loops of size {k,l) to the standard 
Wilson action (one-plaquette action) 






where denotes the average of planar Wilson 

loops of size (/c,Z) and {l,k) in the (p.,^ plane of 
a 4-dimensional lattice of size x . The co¬ 
efficients Ci are, in general, functions of the gauge 
coupling i.e. Ci = Ci{g'^). In the limit g"^ ^ 0 
the two coefficients in Eq. |h should fulfill the rela¬ 
tion Cq^’‘\o) = 1 — (0) in order to insure 

the correct continuum limit. A specific choice of 
^(fc.O(o) allows to eliminate the leading O(a^) cor¬ 
rections on the 

tree level: = -2/[k^l^(k^ + - 2)] 

In the following we will present results for the 
tree level improved actions and which 

have been used recently for thermodynamic calcu¬ 
lations til- These actions can be further im¬ 
proved by adding either additional loops, which 
would allow a systematic perturbative improve¬ 


ment in higher orders of g^, or they may be im¬ 
proved non-perturbatively by decorating the ex¬ 
pansion coefficients with tadpole improvement fac¬ 
tors ll^ 


J /^,0 


( 0 ) 


(k,l) 
’'I,tad 


= U, 


-2(fc+i-2) (fc,i) 


( 0 ) 


Here the factor uq is determined self-consistently, 
for instance from plaquette expectation values , 
^0 = 6l^ Tr □ ). 

Studies at T = 0 have shown that this approach 
indeed leads to a strong reduction of the cut-off 
dependence even on rather coarse lattices. This is, 
of course, most efficient for observables which are 
sensitive to short distance scales. We will discuss 
here in how far this improvement scheme helps to 
improve finite temperature observables. 

Another way of generating an improved action 
within the class of actions defined by Eq. ^ has 
been suggested by Y. Iwasaki [|^ and has been 
extensively used recently for thermodynamic stud¬ 
ies Here renormalization group (RG) meth¬ 

ods have been used to generate an action close to 
the renormalized trajectory, which then has been 
projected onto the two parameter space defined by 
Eq. 1^. This action is defined by the choice 

= cf° = -0.662 . (2) 


Finally we should mention the class of fixed point 
actions, which has been constructed recently 
and which also have been used for thermodynamic 
calculations [|^,^. The fixed point action repro¬ 
duces the continuum dispersion relation and there¬ 
fore does, by construction, also eliminate all cut¬ 
off dependences. These do return, when the fixed 
point action, which depends on an infinite number 
of parameters, is projected on to a limited param¬ 
eter space. A frequently used truncated version is 
expressed in terms of powers of two basic Wilson 
loops, the plaquette, and a non-planar six- 

link Wilson loop, which we denote here by 
The action then reads 

n—l 

k 

+ E E^b" > (3) 

x,u>fi>p n=l 

here Wp,^_p(a:) denotes the average over the four 
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distinct paths of length six on the boundary of a 
3-d cube. In the high temperature, ideal gas limit 
these truncated fixed point actions are indeed very 
similar to the ones defined by Eq. |^. The ideal gas 
limit receives contributions only from the n = 1 
terms in the above sums, ie. it is controlled by 
the two couplings do.i and dip only. The ansatz 
given by Eq. ^ may also be used to define a tree 
level improved action by choosing k = 1 and dop = 
1/3, dip = 1/3. This is similar to Eq. with the 
planar six-link loop replaced by a twisted six-link 
loop. 


3. The High Temperature Ideal Gas Limit 

The importance of improved actions for thermo¬ 
dynamic calculations becomes evident immediately 
from an analysis of the high (infinite) temperature 
limit for a gluon gas on lattices with temporal ex¬ 
tent Nt. In the standard Wilson formulation the 
0{a'^) cut-off dependence in physical observables is 
reflected in large 0{N~‘^) deviations from the con¬ 
tinuum Stefan-Boltzmann limit. On lattices with 
temporal extent Nr one finds, for instance, for the 
deviation of the energy density, €{Nr), from the 
continuum result, esB, ||]0 


<Nr) 

esB 



The large 0{a^) corrections are eliminated by us¬ 
ing the improved actions defined in Eq. ^ One 
indeed finds a strong reduction of the cut-off de¬ 
pendence relative to the standard Wilson formula¬ 
tion. In fact, although one has no control over the 
sub-leading corrections when setting up the sys¬ 
tematic 0{a^) improvement it turns out that in 
many cases not only the 0{a^) corrections get elim¬ 
inated but also the 0{a‘^ corrections get reduced. 
Eor instance, one finds ||] 


improvement does not affect the thermodynamics. 
In how far it leads to an improvement at finite tem¬ 
peratures will be discussed in the next section. 




Figure I. Cut-off dependence of the energy density 
and pressure in the infinite temperature, ideal gas 
limit. The upper figure shows results for some ac¬ 
tions as function of the temporal lattice extent Nr ■ 
In the lower figure results are shown for Nr = 4. 
For further discussions see the text. 


= l + (5) 

with Cl = 0.044(2) for the action and c/ = 

0.178(2) for This leads to a drastic reduc¬ 

tion of cut-off errors as can be seen from Fig. |^. Of 
course, in this infinite temperature limit tadpole 

*’The 0{Nr‘^) coefficient has been given incorrectly in |^. 
It should read 2/5 rather than 2/6. We thank T. Scheideler 
for pointing this out to us. 


In the general case of actions, which also include 
non-planar loops, it is more suitable to calculate 
the pressure, p, (or free energy density, / = —p) 
rather than the energy density. The former is pro¬ 
portional to the logarithm of the partition function, 
and thus is easily obtained in perturbative calcu¬ 
lations. In the lower part of Fig. we also show 
the ratio p{Nr)/psB for Nr = 4 for several actions 
defined through Eqs. and with varying coeffi¬ 
cients Cl and diq, respectively. This also includes 
the case of the RG-improved action defined through 
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Eq. 1^ (open square). As can be seen this action 
leads to rather large cut-off dependent corrections 
in the high temperature limit. In this figure we 
also show results for two versions of the fixed point 
action^ (triangles) and the tree level improved 
actions including planar (open circle) and twisted 
(filled circle) six-link loops. As can be seen all these 
actions do lead to small deviations from the contin¬ 
uum result already on lattices with temporal extent 
Nr =4. 

A similar systematic improvement as discussed 
above for the gluonic sector can be achieved in the 
fermion sector [|^ . We will present some results for 
four-flavour QCD in Section 5. In the next section 
we start with a discussion of results obtained for 
the thermodynamics of the pure SU(3) gauge the¬ 
ory using mainly tree level and tadpole improved 
actions M- 

4. SU(3) Thermodynamics 

Only recently calculations with the standard 
Wilson action could be extended to lattices with 
sufficiently large temporal extent {Nr = 6 and 8) 
that would allow an extrapolation of lattice results 
for bulk thermodynamic quantities to the contin¬ 
uum limit 1^,1^. Computationally the step from 
Nr = 4 to Nr = 8 is quite non trivial as the 
computer time needed to achieve numerical results 
with the same statistical significance on a two times 
larger lattice increases roughly like It therefore 
is highly desirable to use improved actions, which 
suffer less from discretization errors, also for ther¬ 
modynamic calculations. 

We have analyzed the thermodynamics of the 
SU(3) gauge theory using the 0{a‘^) tree level 
(mq = 1) and tadpole {uq < 1) improved actions 
defined in Eq. |^. As discussed above it seems to be 
plausible that improved actions will help to reduce 
the cut-off dependence of thermodynamic quanti¬ 
ties in the high temperature plasma phase. It is, 
however, less evident that this also is of advantage 
for calculations close to the deconfinement phase 
transition where the physics is strongly dependent 
on contributions from infrared modes. In addition 
to a calculation of the equation of state an analysis 
of the cut-off dependence of the critical tempera¬ 
ture itself as well as properties of the first order de¬ 
confinement transition like the surface tension and 

use the notation FP I and FP III for the two sets of 
fixed point actions defined in ||l^. 


latent heat are therefore of interest. 

4.1. The critical temperature 

When comparing the cut-off dependence of ther¬ 
modynamic observables one, of course, has to make 
sure that these are compared at the same value of 
the cut-off. On lattices with fixed temporal extent 
this amounts to a comparison at the same value of 
the temperature, which, for instance, can be de¬ 
fined in terms of the critical temperature for the 
deconfinement transition, i.e. T/Tc- Also the de¬ 
termination of this temperature scale is influenced 
by the finite lattice spacing and will contribute to 
the overall cut-off dependence of thermodynamic 
observables. An indication for the size of cut-off 
dependence in the definition of a temperature scale 
can be deduced from a calculation of the critical 
temperature in units of t/a. The ratio ^fa has 
been studied in quite some detail at different values 
of the cut-off for the SU (3) gauge theory with the 
Wilson action. The 0{a^) dependence has clearly 
been seen in these calculations. An extrapolation 
of these results to the continuum limit yields^ 

^ = 0.631 ± 0.002 . (6) 

Va 

The reduction of the cut-off dependence may best 
be discussed by comparing calculations of T^j 
at a fixed value of the cut-off, for instance for 
oTc = 1/4, with the above continuum extrapola¬ 
tion. Such a comparison is performed in Table 
for some tree level and tadpole improved actions as 
well as the RG-improved action. While in the case 
of the Wilson action Tdy/a deviates by roughly 
10% from the continuum extrapolation at this value 
of the cut-off, the tree level and tadpole improved 
actions only show little deviations. We also note 
that the result for the tadpole improved action is 
consistent within errors with the tree level result. 
The RG-improved action does lead to a somewhat 
larger value for Tdda than all the other calcula¬ 
tions presented in Table |^. However, when com¬ 
paring calculations at different values of the cut-off 
(Fig. g) we note that the RG-improved action as 
well as the other improved actions shows no signif¬ 
icant cut-off dependence, unlike the Wilson action 

“^This differs slightly from the value published in [p^. The 
difference results from a re-analysis of the string tension 
at /3c(A^r) for -/Vr = 8 and 12 which now is based on 
an interpolation of results obtained at nearby values of the 
gauge coupling rather than a single calculation performed 
at l3c{Nr) as it has been done in ]l9|. 
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Table 1 

Critical temperature in units of \f5 on lattices with temporal extent = 4, be. at a value of the cut-off 
given by aT^ = 0.25. Infinite volume extrapolations for the critical couplings have been performed in all 
cases. Further details on the data can be found in the references given in Figure 2. 


action 

/3c 

Tc/^/a 

standard Wilson 

5.69254 (24) 

0.5983 (30) 

(2,2) (tree level improved) 

4.3999 (3) 

0.625 (4) 

(1,2) (tree level improved) 

4.0730 (3) 

0.633 (3) 

(1,2) (tadpole improved) 

4.3525 (5) 

0.634 (3) 

(1,2) (RG improved) 

2.2879 (11) 

0.653 (6)(1) 


where the Oia}) cut-off dependence is clearly visi¬ 
ble. It thus seems that the discrepancy of about 4% 
between the different improved action calculations 
is mainly due to differences in the determination 
of the string tension from the heavy quark poten¬ 
tial and not due to different 0{a^) corrections, be. 
the accuracy is at present limited through ambi¬ 
guities in extracting the long-distance physics from 
the heavy quark potential. This ambiguity is re¬ 
flected in the spread of current determinations of 
the critical temperature, 

Tc/v^ = 0.631 (2) [|§|2|] - 0.656 (4) 0. 



Figure 2. The critical temperature in units of the 
square root of the string tension for various actions 
versus the square of the cut-off. For the tree level 
improved (l,2)-action results for the critical cou¬ 
pling on finite spatial lattices with temporal extent 
Nt = 5 and 6 have been taken from Ref 11 and 
extrapolated to infinite volume fl^ . 


4.2. Bulk thermodynamics 


Bulk thermodynamic quantities like the energy 
density (e) or pressure (p) calculated with the Wil¬ 
son action on lattices of size N^Nr with Nr = 4, 6 
and 8 [0 show a strong cut-off dependence in the 
plasma phase. A first analysis with the tree level 
improved (2,2)-action has shown that this cut¬ 
off dependence gets drastically reduced already on 
an Nr = 4 lattice. This is true also for the tree level 
improved (I,2)-action, the tadpole improved (1,2)- 
action |n| and a fixed point action [0 which all 
yield results consistent with the continuum extrap¬ 
olation performed for the standard Wilson action. 
As an example we discuss here the pressure. 


The pressure can be obtained from an integration 
of the difference of action densities at zero (S'o) and 
finite (St) temperature 0, 


p P 
^ Po 


Nt 



d/3'(.So - St) , 


( 7 ) 


where the zero temperature calculations are per¬ 
formed on a large lattice of size N'^ and the finite 
temperature calculations are performed on lattices 
of size N^Nr- In order to compare calculations 
performed with different actions one has to deter¬ 
mine a physical temperature scale. This can, for 
instance, be achieved through a calculation of the 
string tension at several values of the gauge cou¬ 
pling (/3) and at the critical coupling for the phase 
transition on a lattice of temporal size Nr- This 
yields T/Tc = yPa{(3c)/^/a{j3). 

A systematic analysis of the cut-off dependence 
of bulk thermodynamic observables has been per¬ 
formed previously for the Wilson action |^. Cal¬ 
culations on lattices with temporal extent Nr = 
4, 6 and 8 could be used there to perform an extrap¬ 
olation to the continuum limit. In Fig. ^ we show 
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the results of a calculation of the pressure using 
various improved actions on N^- = 4 lattices. The 
strong cut-off dependence for the Wilson action is 
clearly seen in the upper plot. While all calcula¬ 
tions with improved actions are close to the contin¬ 
uum extrapolation obtained from the Wilson action 
(solid line), the Wilson action result obtained di¬ 
rectly on an Nr = 4 lattice shows strong deviations 
(upper curve). In fact, it seems that the scattering 
of the various improved action calculations around 
the continuum extrapolation mainly reflects the un¬ 
certainties in fixing the temperature scale through 
calculations of the string tension rather than sys¬ 
tematic differences in the cut-off dependence. This 
is also supported by the lower part of Fig. || where 
we show a comparison of improved action calcula¬ 
tions at different values of the cut-off. Calculations 
on Nr = 3 and 4 lattices do not show any signif¬ 
icant cut-off dependence for the tree level as well 
as for the tadpole improved actions. On the other 
hand, the spread of results obtained with different 
actions is larger than the cut-off dependence visible 
for a given improved action. The results also agree 
well with a calculation of the pressure performed 
with a fixed point action on lattices of size 12^ x 3 
(triangles) [|^. It thus seems that systematic er¬ 
rors in the calculation of the equation of state which 
result from the lattice discretization are well under 
control for the SU(3) gauge theory. 

In Fig. ^ we show the energy density, pressure 
and entropy density obtained from an extrapola¬ 
tion of results obtained with the standard Wilson 
action. We note the rapid rise in all observables at 
Tc followed by a rather slow approach to the asymp¬ 
totic ideal gas behaviour. The latter is in accor¬ 
dance with the expectation that the high tempera¬ 
ture behaviour of QCD is controlled by a universal 
function, which only depends on a running coupling 
that varies logarithmically with temperature. How¬ 
ever, although e, s and p rapidly come close to the 
ideal gas limit, the (15-20)% deviations observed at 
temperatures as large as 5Tc are still too large to 
be described by perturbation theory. The pertur¬ 
bative expansion of the thermodynamic potential 
is converging quite badly and would require a run¬ 
ning coupling which is significantly smaller than 
unity [^3y24| ] . As a consequence, the perturbative 
expansion seems to converge only for temperatures 
much larger than Tc. In particular, the calculation 
of screening lengths on the lattice do, however, sug¬ 
gest that the running coupling is larger than unity 
even at 5Tc [p^ . 




Figure 3. Pressure of the SU(3) gauge theory cal¬ 
culated with the Wilson action and different im¬ 
proved actions on Nr = 4 lattices (upper figure). 
The lower figure shows a comparison of calcula¬ 
tions with tree level and tadpole improved actions 
on Nr = 3 and 4 lattices. Also shown there are 
results from a calculation with a fixed point action 
(triangles). The arrows indicate the ideal gas result 
in the continuum limit. 


4.3. Surface tension and latent heat 

The success of improved actions for the calcula¬ 
tion of bulk thermodynamics even at temperatures 
close to Tc naturally leads to the question whether 
these actions also do lead to an improvement at Tc 
itself where a strong cut-off dependence has been 
observed previously in the calculation of the latent 
heat (Ae) and the surface tension (cr/) in studies 
with the Wilson action on lattices with temporal 
extent up to = 6 [|^J2^. The region around 
Tc is, of course, a highly non-perturbative regime. 
However, observables like Ae and cr/, which char¬ 
acterize the discontinuities at the first order decon¬ 
finement phase transition in a SU (3) gauge theory. 
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Figure 4. Extrapolation to the continuum limit for 
the energy density, entropy density and pressure 
versus T/T^- The dashed horizontal line shows the 
ideal gas limit. The hatched vertical band indicates 
the size of the discontinuity in e/T"' (latent heat) 
at Tc- Typical error bars are shown for all curves. 


do depend on properties of the low as well as the 
high temperature phase. As the latter is largely 
controlled by high momentum modes it may be 
expected that some improvement does result even 
from tree level improved actions. 


We have extracted cr/ |n| from the probability 
distribution of the absolute value of the Polyakov 
loop, P{\L\), following the analysis presented in 
Ref. 1^^. The probability distribution at the min¬ 
imum is proportional to 


P(|L|) - exp(- [/iPi + / 2 P 2 + 2aiA] /T) (8) 


where fi denotes the free energy in the phase *, 
Vi is the volume occupied by that phase and A 
denotes the interface area of the finite system. The 
distribution functions for Polyakov loop at Tc show 
a double peak structure. From the depth of the 
minimum between these peaks one can extract the 
surface tension. 

In Fig. P we show results for the surface tension 
calculated with the Wilson action as well as tree 
level and tadpole improved (l,2)-actions at differ¬ 
ent values of the cut-off. It is quite remarkable that 
the results for the tadpole improved action show 
practically no cut-off dependence when comparing 
calculations on = 3 and 4 lattices. Furthermore 
it is evident that both improved actions lead to sig¬ 


Figure 5. Comparison between surface tension cal¬ 
culations performed with the Wilson action as well 
as tree level and tadpole improved (l,2)-actions. 


nificantly smaller values for the surface tension at 
finite values of the cut-off. In fact, the results for 
Nt = 4 are in agreement with a quadratic extrap¬ 
olation of the Wilson data to the continuum limit. 

In Table || we give results for aj on the largest 
lattices considered. Clearly the surface tension ex¬ 
tracted from simulations with improved actions on 
lattices with temporal extent N^. = 4 are substan¬ 
tially smaller than corresponding results for the 
Wilson action. In fact, they are compatible with 
the Nr = 6 results for the Wilson action. While the 
cut-off dependence of the surface tension calculated 
on finite lattices seems to be well under control now 
through the use of improved actions, it remains to 
get control over the behaviour in the thermody¬ 
namic limit 00 ). The difference between 

the surface tension calculated on the largest spatial 
lattice and the extrapolated value is about 20%. 
However, it is not quite clear whether the ansatz 
used to perform the extrapolation to the infinite 
volume is already applicable for the currently used 
spatial lattice sizes. It includes, for instance, the 
contribution of translational zero modes, which are 
not yet observed in the present calculations j^. 
Given these caveats we find from an extrapolation 
of the results obtained with the tadpole improved 
actions [|ll[ 

^ = 0.0155 ±0.0016 (9) 

^ C 


The latent heat is calculated from the disconti- 













Table 2 

Surface tension and latent heat for improved actions and the Wilson action. Results for the Wilson action 
are based on data from using the non-perturbative /3-function calculated in [|^ . Improved action results 
are taken from ||ll|]. 


action 


Nr 


Mn 

standard Wilson 

24^ X 36 

4 

0.0300 (16) 

2.27 (5) 


302 X 48 

6 

0.0164 (26) 

1.53 (4) 

(l,2)-action (tree level) 

32^^ 

4 

0.0116 (23) 

1.57 (12) 

(I,2)-action (tadpole) 

323 

4 

0.0125 (17) 

1.40 (9) 


nuity in (e — 3p). This in turn is obtained from the 
discontinuity in the various Wilson loops entering 
the definition of the improved actions. In the case 
of actions, which are defined through a set of cou¬ 
plings, which do not depend on the gauge coupling 
g'^, the latent heat is, in fact, simply proportional 
to the discontinuity in the action expectation value 
at Tc- In general one has, however, to take into 
account contributions resulting from derivatives of 
the couplings with respect to The latent 

heat is then given by 


Ae 



1 

6 





( 10 ) 


with S = S — dS/df3. 

The discontinuity in the expectation value of S at 
Pc is obtained by calculating these separately in the 
two coexisting phases at Pc- In order to extract the 
latent heat one still needs the /3-function entering 
the definition of Ae/T^ in Eq. |^. The necessary 
relation a(/3) has been obtained from a calculation 
of ^/aa (improved actions) or a determination of 
TcQ (Wilson action) |^. Results for on 

Nr = 4 and 6 lattices are summarized in Table 

A comparison with simulations on even coarser 
lattices (W = 3) shows that the cut-off depen¬ 
dence in the latent heat follows a similar pattern 
as that shown for the surface tension in Fig. H. In 
particular, the results obtained with the tadpole 
improved action on Nr = 3 and 4 lattices agree 
with each other within errors. We thus may use 
these results to estimate the latent heat in the con¬ 
tinuum limit. Using Eq. and = 420 Mev we 
find 


Ae ~ I.5T^^ = 0.23cr^ = 0.9 GeV/fm^ 


5. Four-flavour QCD with an improved stag¬ 
gered action 

In the fermionic sector of QCD the influence of 
a finite cut-off on bulk thermodynamic observables 
is known to be even larger than in the pure gauge 
sector. For instance, in the staggered formula¬ 
tion the energy density of an ideal fermi gas dif¬ 
fers by more than 70% from the continuum value 
on a lattice with temporal extent Nr = 4 and ap¬ 
proaches the continuum value only very slowly with 
increasing Nr- This cut-off dependence can dras¬ 
tically be reduced with an O(a^) improved stag¬ 
gered action. In a first attempt to analyze the im¬ 
portance of improvement in the fermion sector we 
have performed calculations with an improved ac¬ 
tion, S^[U] = where a higher order 

difference scheme (one-link and three-link terms), 
is used to improve the fermionic part [^ j. The 
improved fermion matrix reads 

M[U]xy = mSx,y + ^ —A [U]^ y — —33 [U]^ 
with 

A[[/]a;y = r]x,y, {Ux,yAx,y-jl — Ul_p^ y^Sx,y+p^ 

B[U] xy — E Vx,y,{Ux^fi Ux-\-fL,fiUx-\-2ft,fi^x,y—3il 

-[/■'' - Sxv+3u) 

where r]x,y. = (—and r]x,o = 1 de¬ 
note the staggered phase factors. In the gluonic 
sector we use the tree level improved (l,2)-action. 
With this action the overall cut-off distortion of 
the ideal gas limit on a 16^ x 4 lattices reduces 
to about 20%. We have performed simulations for 
two quark masses, ma = 0.05 and O.I H- Like 
in the pure gauge sector the improvement is visible 
already close to Tc- 


( 11 ) 
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The general structure of the four-flavour equa¬ 
tion of state does not differ much from that of the 
pure gauge theory. In fact, the temperature depen¬ 
dence of the pressure is very similar to that of the 
pure gauge theory, if we re-scale the latter by an 
appropriate ratio of the number of degrees of free¬ 
dom so that the high temperature limit coincides 
for both cases. This is shown in Fig. p. 



0.5 1 1.5 2 2.5 3 3.5 


Figure 6. Comparison between the pressure of four- 
flavour QCD on a 16^ x 4 lattice for two values of 
the quark mass and the pressure of the pure SU (3) 
gauge theory. The latter has been re-scaled by the 
appropriate number of degrees of freedom of four- 
flavour QCD (29/8) and a factor 1.19 which takes 
care of the remaining cut-off distortion of the ideal 
gas limit resulting from the use of the Naik action. 


While the pressure is obtained in complete anal¬ 
ogy to the pure gauge theory through an integra¬ 
tion of action differences (Eq. |^), the calculation 
of the energy density now also requires the deter¬ 
mination of differences of the chiral condensates at 
zero and non-zero temperature as well as the cut¬ 
off dependence of the two bare couplings, f3 and 

nig, 


e — 3p 

2^4 


-iVi 


d/3 
din a 




dma 


din a 




( 12 ) 


Only in the chiral limit the derivative dlnma/da 
vanishes and (e — 3p) is again proportional only to 
the /3-function, d/3/dlna, as it is the case in the 
pure gauge sector. 


Results for the energy density are shown in 
Fig. 1^. The energy density does stay close to the 
ideal gas limit immediately above T^. We do ob¬ 
serve an overshooting of the ideal gas limit close 
to Tc for the non-zero quark masses considered 
by us. This is a feature not seen before in the 
pure gauge sector. Whether this will persist for 
finite values of the quark mass or is an artifact 
of our present statistical accuracy has to be clar¬ 
ified in further more detailed investigations. The 
overshooting seems, however, to disappear in an 
’’extrapolation” to the chiral limit, which we con¬ 
structed by ignoring the term proportional to (^/'0) 
in the definition of the energy density |^. The 
contribution of this term vanishes in the chiral limit 
as the derivative dlnma/da is proportional to the 
quark mass. 



Figure 7. Energy density of four-flavour QCD on 
a 16^ X 4 lattice. The lower set of curves shows an 
’’extrapolation” to the chiral limit which has been 
obtained by ignoring the second term in Eq. 10 (see 
text). 


We note that the energy density in the critical 
region is larger than in the pure gauge case when 
expressed in units of Tc, i.e. e ~ lOT^. However, 
as the critical temperature is substantially smaller 
than in the pure gauge theory, the critical energy 
density again turns out to be similar, i.e. about 
1 GeV/fm3. 

6. Improved Fermion Actions 

The calculations with an improved fermion ac¬ 
tion show that a strong reduction of the cut-off de¬ 
pendence is possible in the high temperature phase. 
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Table 3 

Coefficients for rotationally invariant staggered 
fermion actions 


{j,k,l,m) 

Cj,k,l,m 


p4-act ion 

p6-action 

(I,0,0,0) 

0.375 

0.32 

(I,2,0,0) 

0.0208333 

0.02 

(I,2,2,0) 


0.0010938 

(3,0,0,0) 


0.0047917 

(1,2,2,2) 


0.00125 

(3,2,0,0) 


0.00125 


Still a further improvement is necessary in order to 
reduce the cut-off dependence to only a few percent 
as it is the case in the pure gauge sector. 

Although the Naik action significantly reduces 
the cut-off dependence of thermodynamic quanti¬ 
ties calculated with the staggered fermion action 
the deviations from the continuum result are still 
about 20% on lattices with temporal extent Nt = 
4. Moreover it seems that the flavour symmetry, 
which is broken at C>(a^) in the staggered fermion 
action, is not much improved in the Naik formu¬ 
lation 11^. It thus is desirable to look for more 
suitable fermion actions. In general, we may dis¬ 
cretize the free fermion action not only by introduc¬ 
ing nearest neighbour couplings between fermion 
fields but we may also allow for couplings between 
sides further apart, 

*S'_F = ^ ^ Cj^k,l,m ■ 

■ [tp{x + ja^ + kai, + lap + maa) 

—ip{x — jup — fcoj/ — lap — macr)] ■ (13) 

A specific choice of couplings Cj^k,i,m has been given 
in for a staggered fixed point action. This does, 
however, include too many non-zero couplings for 
being useful in a numerical simulation. 

From the analysis of the improved gluon actions 
involving six-link operators we learned that non- 
planar loops are, in fact, more efficient in reducing 
the cut-off dependence, even when the 0{a^) cor¬ 
rections are not eliminated exactly. One thus may 
consider to replace the straight three-link term in 
the Naik action by a bended three link term. The 
coefficient of this term can be fixed by demand¬ 
ing that the fermion propagator is rotationally in¬ 
variant up to 0{p‘^). Similarly one can construct 


actions which include longer paths and are rota¬ 
tionally invariant up to 0{p^). We call these ac¬ 
tions p4 and p6-action, respectively. For instance, 
in order to obtain a fermion propagator, which is 
rotationally invariant up to 0(j)^) within the class 
of staggered fermion actions that contain non-zero 
couplings only for one and three-link paths, the co¬ 
efficients are constrained by 

Cl,0,0,0 + 3 c3,0,0,0 + 6ci,2,o,o = 2 

Cl,0,0,0 + 27 c3 0,0,0 + 6ci_2,o,o = 24ci 2,0,0 • (14) 

Similar constraints can be obtained for actions that 
are rotationally invariant up to 0{p^). A possible 
set of non-zero couplings for these actions is given 
in Table 3 and a calculation of the fermion energy 
density is displayed in Fig. ||. We note that with 
the simple choice of the p4-action, which includes 
in addition to the standard one-link term only a 
bended three-link term, the cut-off dependence is 
reduced to a few percent even on Nr = ^ lattices. 
As the introduction of gauge fields in this action 
will require smeared paths like in the fat link action 
proposed in there is the possibility that such 
actions also improve the flavour symmetry of the 
staggered fermion formulation. This, however, has 
to be examined still in more detail. 



Figure 8. Cut-off dependence of the fermion energy 
density calculated from several improved actions. 


7. Conclusions 

Thermodynamic observables of the SU{3) gauge 
theory and QCD studied with improved gauge and 
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fermion actions show a drastic reduction of the cut¬ 
off dependence in the high temperature limit as well 
as at Tc- The major improvement effect is already 
obtained with tree level improved actions. The cal¬ 
culation of the equation of state of the SU (3) gauge 
theory seems to be well under control and major 
sources for lattice artifacts have been eliminated. 
At least for finite temperature calculations on lat¬ 
tices with temporal extent Nt = 4 or larger system¬ 
atic effects due to additional tadpole improvement 
are within the current statistical accuracy. Tadpole 
improvement leads to further reduction of system¬ 
atic errrors on N^. = 3 lattices. This opens the 
possibility to repeat quantitative studies of various 
other thermodynamic quantities on rather coarse 
lattices which so far could only be investigated on 
a qualitative level without reliable extrapolations 
to the continuum limit. 

The calculation of thermodynamic quantities 
with improved fermion actions has just started. In 
this case further work is still needed in order to be 
able to select an appropriate action which even in 
the ideal gas limit has as little cut-off dependence 
as the tree level improved gluon actions. 
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